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Abstract 



ing states in the He system employing a realistic nucleon-nucleon potential in 
the framework of the resonating group model (RGM). We present for comparison 

3 : 

with these microscopic RGM calculations the results from a charge- independent, 



Coulomb-corrected i?-matrix analysis of all types of data for reactions in the A = 4 
system. Comparisons are made between the phase shifts, and with a selection of 
measurements from each reaction, as well as between the resonance spectra obtained 
from both calculations. In general, the comparisons are favorable, but distinct dif- 
ferences are observed between the RGM calculations and some of the polarisation 
data. The partial-wave decomposition of the experimental data produced by the R- 
matrix analysis shows that these differences can be attributed to just a few 5-matrix 
elements, for which inadequate tensor-force strength in the N — N interaction used 
appears to be responsible. 

PACS numbers: 21.40. +d, 21.60.-n, 25.10.+S 
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1 Introduction 



The 4 He atomic nucleus is one of the best studied nuclei, both experimentally and theo- 
retically, as summarized in the recent A = 4 compilation IJ. Besides the many textbook 
examples of the gross structure, there are subtle points, usually neglected, which yield 
large effects. Most of these effects are only qualitatively and in most cases never quan- 
titatively explained. None of the existing calculations aims at a complete understanding 
of the many features of 4 He no wonder in view of the many different phenomena studied 
so far rjTJI . With the recent compilation |TJ and the comprehensive .R-matrix analysis 
of a large amount of scattering data below 10 MeV, a microscopic calculation for the 4 He 
system in this energy range is most needed for. 

The difficulties arise already for the total energy of the ground state. In most cases, 
parameters of the effective nucleon-nucleon potential are modified in such a way as to 
reproduce that energy or the threshold energies, but seldom both simultaneously. Calcu- 
lations using effective potentials, however, allow to reproduce only some features, while 
failing in others. Therefore it remains unclear if the failure is due to the deficiencies of 
the potential or of the model itself. Only calculations employing realistic nucleon-nucleon 
potentials without any adjustable parameters can help to answer this question. Such 
calculations, however, are very complicated and very time consuming, no matter what 
type of approach is chosen. For realistic potentials, calculating the ground state energy is 
already a major effort and only three-body forces allow to reproduce that energy within 
a few keV in Greens-function Monte-Carlo calculations [[J. The next open question is the 
structure of the first excited state, the + state just between the 3 H — p and the 3 He — n 
thresholds, considered frequently to be a breathing mode, which is reasonably narrow. In 
calculations this state quite often turns out to be bound (see e.g. [|]). On the experimen- 
tal side this state shows a peculiar behaviour, indicating a very broad structure in a — a 
breakup ||. The i?-matrix analysis |2j indicates the possible origin of this feature. All 
the other resonances do not exihibit any narrow structure, but turn out to have a decay 
width of the order of MeV [|T], ^| , which makes their location rather arbitrary. 

On the other side there is such an amount of scattering data, that phase shift analyses are 
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possible, especially for proton-triton and 3 He-neutron scattering. The i?-matrix analysis 
[[J connects all the possible reaction channels simultaneously, thus giving the most com- 
plete description and allowing to interpolate to every desired energy or reaction. Since 
the number of parameters in this analysis is in the order of one hundred, however, an 
unrestricted fitting procedure is not possible and some physical input is needed, which 
will be discussed below. 

The aim of this paper is to show that with present-day computers, a rather unrestricted 
calculation using realistic nucleon-nucleon potentials, like the Bonn potential Jj]], is feasible 
and yields essentially all known experimental features. What we consider most surprising 
is the close agreement between the i?-matrix analysis, which has as input only data in the 
4 He system, and the refined Resonating Group calculation ||, which has as input only a 
realistic nucleon-nucleon potential without any adjustable parameters. 

We organize the paper in the following way: the next chapter contains a brief review 
of the Resonating Group Model (RGM) and the model space used for the calculation, 
where we want to point out the essential differences with existing calculations. After 
that, we review the R- matrix analysis with special emphasis on all restrictions employed 
in the fitting of the data. The next chapter is devoted to a detailed description of the 
calculation for the ground state of 4 He. Then we start with a comparison of the RGM 
and i?-matrix results for the various partial waves. Finally we compare selected data with 
the results of the i?-matrix analysis and the RGM-calculation, and give a brief outlook. 



2 Resonating Group Model and the Model Space for 
4 He 

The Resonating Group Model in its various modifications || ^, |H| is a suitable method to 



calculate the scattering of composite objects. In this case the main technical problem lies 
in the evaluation of the many-body matrix elements. In order to facilitate the calculation 
of the orbital matrix elements, however, all radial dependencies have to be of Gaussian 
form, since only then proper antisymmetrization for the translationally invariant wave 
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function is feasible. Therefore also the potential has to be given in terms of Gaussian 
functions. In this work we use the r-space version of the Bonn-potential |7[] expressed in 



Gaussian functions | IT| . All further reference to the Bonn-potential is indeed a reference 



to this Gaussian version. 

Powerful techniques have been developed to achieve the analytic calculation of the indi- 
vidual matrix elements 0, [K| . Even with these techniques, calculations are usually far 



from trivial, hence most calculations are restricted to single channel problems, neglecting 
e.g. effects of the tensor-force, and/or use rather simple effective potentials. Since the 
calculation of the orbital matrix elements is the most time consuming part, symmetries 
of the orbital wave function are exploited as far as possible M . For simple (central or 
non-central) potentials, single harmonic oscillator wave functions yield already a good 
description of the ground state wave function of the lightest nuclei (see e.g. [|J). For a 
realistic N — N potential this is no longer the case, because of the strongly repulsive short- 
ranged core of the potential, which suppresses the wave function for small internucleon 
distances. In many cases even a node occurs in the wave function at small distances. In 
terms of Gaussians this means that more than one Gaussian width parameter is needed 
to describe even the simplest wave function. Furthermore, the tensor force plays a crucial 
role. 

Let us explain the situation for the simplest nucleus, the deuteron: For an effective 
N — N interaction a single S'-wave binds the deuteron already; for a reasonable radius of 
the wave function a linear combination of two Gaussian functions is enough (see [||, [5|]). 
For a realistic potential only a superposition of S- and D-waves yields a bound deuteron, 
with at least 3 Gaussians on the S'-wave and two on the D-wave for the Bonn potential 
[0] (see Appendix A for the detailed wave function). For other potentials, e.g. [O, we 
find the same behaviour. 

In the 4 He system we have three two-fragment channels, the triton-proton, the 3 He- 
neutron, and the deuteron-deuteron channel. For the latter, the D-w&ve admixture in 
the wave function leads to an additional coupling of internal angular momenta of 2 to 
the relative orbital angular momentum. These couplings lead to very complicated wave 



3 



functions for the composite system and increase the computing time neccessary by orders 
of magnitude. Since we allow for orbital angular momenta inside the deuteron, we treat 
the 4 He system essentially as a four nucleon system, i.e. we consider it as a four cluster 
problem in the framework of the RGM. 

Also the triton-proton channel is treated as a four-cluster problem, because we need 
internal D-waves to bind the triton. To get a reasonable binding energy, we need D- 
waves on all internal coordinates, thus leading to a three-cluster description of the triton. 
Here we need three different Gaussian width parameters on both internal coordinates. To 
lower the binding energy even more, we allow for a second set of width parameters for 
all components of the wave function containing at least one D-w&ve. Some details are 
given in Appendix A. Since in 3 He only a neutron is exchanged against a proton compared 
to the triton, the structure of the wave function is identical, with just modified sets of 
Gaussian width parameters, given in Appendix A. 




Figure 1: Angular momentum structure of the 3 He wave function. 

Allowing for all possible combinations of li < 2 and I2 < 2 (see fig. [I]) that can contribute to 
J n = l/2 + , we get a total of 23 combinations. The total binding energy and the individal 
contributions of kinetic and potential energy compare favourably well with the results 
of Carlson [^] (see table [1]). This wave function, however, is too complicated to be used 
in a full-fledged scattering calculation. Therefore we restricted the triton wave function 
to the dominant terms. These are all three spin-isospin combinations with l\ — I2 — 
and furthermore all combinations with the proton-neutron subsystem having spin=l and 
pairs h, l 2 coupled to l 3 are [1, 1]2, [2, 0]2, [0, 2]2, [2, 2]0, [2, 2]1, and [2, 2]2. This restricted 
wave function consists of only 10 combinations and yields only about 70 keV less binding 
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Table 1: Comparison of triton results with other calculations and 3 He results 



calculation 


energy 


(V) 


<-£a> 




Pltot=2 




[MeV] 


[MeV] 


[MeV] 


[fm] 


[%} 


Carlson (var) 


-7.96(03) 


-48.2(05) 


40.3(05) 


1.69(02) 


7.0(1) 


Carlson(Faddeev) 


-8.29 


-49.0 


40.7 


1.69 


7.0 


full wave fuction 


-8.015 


-48.875 


40.860 


1.64 


6.95 


reduced w.f. 


-7.951 


-48.587 


40.637 


1.64 


6.86 


3 He full w.f. 


-7.343 


-47.464 


40.121 


1.69 


6.89 


3 He reduced w.f. 


-7.258 


-47.124 


39.866 


1.68 


6.76 



energy than the full wave function (see table [I]), but saves large amounts of computing 
time. In test cases we convinced ourselves that this reduction of the wave function of 
our physical channels has only minor effects on the 4 He results. Also for 3 He we use the 
similarly restricted wave function. 

Together with the deuteron wave function we find total binding energies and threshold 
energies which are given in table [| Especially the relative energies compare favourably 
well with the experimental numbers. 

Table 2: Comparison of experimental and calculated total binding energies and relative 
thresholds (in MeV) for the reduced wave functions 



channel 




Ethres 




exp. 


cal. 


exp. 


cal. 


3 H-J9 


-8.481 


-7.951 






3 He — n 


-7.718 


-7.258 


0.763 


0.693 


d — d 


-4.448 


-3.809 


4.033 


4.142 



All further physical channels are three- or four-body breakup channels, which cannot in 



principle be treated within the RGM [13]. To allow for the neccessary amount of flux 
that goes into these channels, we approximate these breakup channels by two-fragment 
channels formed by appropriate combinations of the singlet deuteron d with itself or with 
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a deuteron. The approximate bound state wave function for the singlet deuteron is also 
given in Appendix A. 

This model space, however, is by no means sufficient to find reasonable results in all partial 
waves. We have to add so called distortion channels ||13|| , sometimes called pseudo-inelastic 
channels [fLCR , to allow for enough variational freedom in the interaction region. These 
distortion channels contain therefore no asymptotic part, so that they cannot carry any 
flux, but have e.g. different symmetries or just increase the model space such that certain 
parts of the interaction can yield larger contributions. Some details about the distortion 
channels will be given during the discussion of the individual partial waves. 

Having specified the internal wave functions of our fragments, we indicate the structure 
of the total wave functions; details of the construction of the wave functions can be found 
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TBI. To keep the notation transparent, we give only the gross structure of the wave 



function. 
¥ = A 



[^ 3 H0p] 5 Xm p Oh p ) + [^He^n] 5 xfifen ( n Hen) 
L,S L,S 

+ [^d0d] S Xdd" ( r dd) + distortion channels 

L,S 



(1) 



Here A denotes the four particle antisymmetrizer. The relative-motion wave functions 
are of the structure 



X L (r) = (b L - f L + a L -~g L + Y. c ^ " E »A ' y l(P), 



(2) 



where the fi and gi are regular and regularized irregular Coulomb functions for the 
appropriate channel and the E u l are Gaussian functions of width j3 v times r L , which take 
care of the wave function being different in the interaction region from the asymptotic 
form. Usually the coefficients hi are determined by the appropriate boundary condition, 
i.e. only one element in eq.(l) being unity and all others zero. The coefficients cll and c u l 
are linear variational parameters determined from the Kohn-Hulthen variational principle 



1J via 

( 5^ \H - E\ = 0. (3) 
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S denotes the channel spin and L the relative angular momentum. The S'-matrix is 
directly related to the reactance matrix elements a, via the Caley transformation, S = 
(1 + za) • (1 — ia)^ 1 ; for further details of the method see || [T3|. 



The i?-matrix formulation |18j of multichannel nuclear reaction theory gives a powerful 
and convenient energy- dependent framework for describing the experimental measure- 
ments. It parametrizes the unitary S'-matrix in terms of (real) reduced-width amplitudes 
7 c a and eigenenergies E\, for fixed values of boundary- condition numbers B c and chan- 
nel radii a c . Since these parameters reflect the nature of the interactions only at small 
distances (r c < a c ), they can be approximately constrained by the symmetry properties 
of the strong (nuclear) forces. In this analysis of reactions in the A = 4 system, the 
approximate charge independence of nuclear forces was used to relate the parameters in 
charge-conjugate channels, while allowing simple corrections for the internal Coulomb 
effects. 

The Coulomb-corrected, charge- independent R- matrix analysis of |L5| (also used in 0) 
takes its isospin T = 1 parameters from an analysis of p— 3 He scattering data |16|] that gives 
a good description of all data at proton energies below 20 MeV. The eigenenergies are, 
however, shifted by the internal Coulomb energy difference AEc = —0.64 MeV and the 
p — 3 H and n — 3 He reduced width amplitudes are reduced by the isospin Clebsch-Gordan 
coefficient l/y/2. The isospin T = parameters are then varied to fit the experimental 
data for reactions among the two-fragment channels p — 3 H, n — 3 He, and d — 2 H, at 
energies corresponding to excitations in 4 He below 29 MeV. A summary of the channel 
configuration and data included for each reaction is given in table In this fit, the T = 
nucleon-trinucleon reduced-width amplitudes are constrained by the isospin relation 



and a small amount of internal Coulomb isospin mixing is introduced by allowing 



3 The R- matrix analysis 



7n3He 




(4) 




(5) 
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Table 3: Channel configuration (top) and data summary (bottom) for each reaction in 
the 4 He system i?-matrix analysis 



Channel 


^max 


a c (fm) 


3 H-J9 


3 


4.9 


3 He— n 


3 


4.9 


2 R-d 


3 


7.0 



Reaction 


Energy range 


(MeV) 


# Observable types 


# Data points 


3 H(p,p) 3 H 


E p = 0- 


11 


3 


1382 


3 H(p, n) 3 He + inv. 


E p = 0- 


11 


5 


726 


3 He(n, n) 3 He 


E n = 0- 


10 


2 


126 


2 R{d,p) 3 E 


E d = 0- 


10 


6 


1382 


2 H(d,n) 3 He 


E d = 0- 


10 


6 


700 


2 E{d,d) 2 E 


E d = 0- 


10 


6 


336 




totals: 


28 


4652 



which is neccessary to reproduce the differences between the two branches of the d — d 
reactions. Note that this i?-matrix analysis is not yet completed nor fully documented. 
However, this analysis represents the most comprehensive and detailed attempt to date 
to give a unified phenomenological description of the reactions in 4 He. 

Since the level information is derived from a multilevel i?-matrix parametrization, a gen- 
eralization of the single-level Breit-Wigner prescription is neccessary. For the convenience 
of the reader, we repeat here the essential steps given in . 

The multilevel generalization of the usual single-level prescription in terms of i?-matrix 
parameters is to find the poles and residues of the channel-space matrix 

K R = G[e{E)-E}- 1 G T 



in which the elements of G are 

G cX = P c 1/2 7c a 



(6) 



(7) 
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and e(E) is the level-space matrix of elements 

e X'X = E X 5 X , X - £ ll c (S c - B c ) 7cA . (8) 

c 

In eqs.(7) and (8) S c and P c are, respectively, the usual energy-dependent channel shift 
and penetrability functions. As the notation implies, Kr in eq. (6) is similar to Heitler's 
reactance matrix |pL7|| , except that it is not a true asymptotic quantity. It therefore gives 
resonance parameters that depend on the values of a c , so that properties such as level 
spacings are functions of the channel radii. 

Near one of its poles, Kr has the rank one form 

1 PRPR fQ\ 

in terms of the channel-space residue amplitude pR and pole energy E R . The pole energy 
is taken to be the "resonance energy" and the total width is 

r = ^ = E4, (io) 

c 

naturally suggesting that partial widths be defined as 

t c = pL (ii) 

This prescription gives resonance parameters based on the positions and residues of ap- 
parent poles of the S- matrix, as seen from the real energy axis of the physical sheet. 

When the S"-matrix is continued onto the complex energy surface, near one of its poles it 
has the form 

S = i-^- (12) 

where Eq = Er — iY/2 is the complex pole energy and po is the complex residue amplitude. 
A procedure for obtaining Eq and po from R- matrix parameters is given in [|19|]. The 
expectation of the Breit-Wigner approximation is that on any Riemann sheet, Er will 
have the value given in eq. (9), V will be given by eq. (10) and the residue amplitudes 
Po will differ from pR only by unimportant phase factors. For light systems like the 4 He- 
system described in terms of i?-matrix parameters, however, this is often not the case 
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[|P9| , pCfj . As explained in [19|, a parameter characterizing the strength of an S'-matrix 
pole, 

S P = 4? (13) 
in terms of the magnitude of its residue compared to its displacement from the real axis 
can be quite different from unity. 



Furthermore, "shadow" poles [51] associated with a resonance can have different positions 
and residues on different sheets of the Riemann energy surface due to extended unitarity 
of the S'-matrix. We refer to a resonance exhibiting any of these differences with the 
Breit-Wigner expectations as a non-Breit-Wigner resonance. 

In order to compare the RGM-calculation and the i?-matrix analysis also in the level 
parameters, we adopt the same prescription in both cases. In the i?-matrix analysis 
channel radii a p t = a n 3 He = 4.9 fm and add = 7.0 fm are used. In the RGM calculation we 
"locate" therefore a resonance at that energy, where the difference between the calculated 
phase shift and the "background phase shift", determined from the appropriate channel 
radius, passes through ninety degrees. 



4 The 4 He Ground State 

In this chapter we will give some information about the ground state of 4 He. For this 
purpose, we consider different model spaces. The simplest space consists just of the three 
1 So channels 3 H — p, 3 He — n, and d — d. In such a model space we find only half the 
experimental binding-energy of -28.296 MeV (see table Adding the 5 Do d — d channel 
and the d — d channel to simulate the breakup channels, we gain about 2 MeV (see table 
i.e. taking all the physical channels into account, the ground state is only 8 MeV below 
the lowest threshold. To improve the binding energy, we have to add distortion channels. 
First of all we allow for 30 3 H — p distortion channels and 30 3 He — n ones, which increases 
the binding energy by about 8 MeV (see table [|). Adding 82 deuteron-deuteron distortion 
channels we gain again about 1.5 MeV. Finally we allowed for additional 62 distortion 
channels of the structure 2N — 2N, which brought the energy down to -25.849 MeV, close 
to the result of ref. [[|. 
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Table 4: Calculated ground state energies in MeV for various model spaces 

ref. 1 
-25.86(15) 
81.62(7) 
0.74(01) 



# channels 


3 


5 


5 + 60 


5+60+82 


5+204 


227 


total energy 


-14.033 


-15.873 


-24.345 


-25.798 


-25.849 


-25.910 


kinetic energy 


36.671 


44.834 


74.767 


80.642 


80.834 


80.981 


Coulomb 


0.611 


0.641 


0.772 


0.795 


0.796 


0.797 


Central 


-21.238 


-23.489 


-37.320 


-39.714 


-39.756 


-39.828 


Tensor 


-20.833 


-27.881 


-46.501 


-50.496 


-50.684 


-50.772 


Spin-orbit 


0.115 


0.236 


0.418 


0.474 


0.481 


0.449 


p 2 -potential 


-9.359 


-10.214 


-16.482 


-17.499 


-17.519 


-17.537 



1.04(01) 
-17.77(27) 



Increasing the model space even further improved the binding energy only marginally, 
but needed much more computing time and caused almost numerical linear dependencies; 
therefore we did not pursue it any further. Since e.g. the triton wave function consists 
already of 10 different components, we finally diagonalized the Hamiltonian in the full 
space spanned by all structures calculated and no longer coupled the various components 
to physical channels. The 5 physical channels are spanned by 36 different components. 
We omitted the small width parameters in these components, which are neccessary to 
reproduce the oscillatory behaviour of the scattering wave functions, but correspond to 
large distances between the fragments. Hence they lead to minimal interaction between 
the fragments and thus yield eigenvalues of the Hamiltonian close to the threshold energies. 
We convinced ourselves that a reasonable choice is to choose the same width parameters 
as for the distortion channels. The actual values are given in Appendix B. Because 
of numerical linear dependencies, only 18 out of the 36 components could be retained, 
yielding a total of 227 channels. The results are displayed in the next-to-last column 
of table |4j. The various results are in close agreement with those of ref. [[| . The point 
Coulomb form factor is also quite similar to that of ref. [2~I] found for the Argonne 
potential (see fig. |2|). Therefore we believe that the inclusion of meson-exchange-currents 
into our calculation would also reproduce these data. 

Our calculation with 227 channels corresponds to a calculation in a certain configuration 
space with only bound states. Hence, we consider the two lowest eigenvalues of the 
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10 u p- 



10"' =r 



10^ =- 



10" 



10" 




0.0 



j L 



2.0 



J L 



4.0 



6.0 



q[fm ] 



Figure 2: Point Coulomb form factor for the 4 He ground state (full line) together with 
data [23, [26| and the results of ref. (dashed line). Also shown is the result for the 
minimal wave function (dotted line). 



Hamiltonian in this space as variational approximations to the ground and first excited 
+ states. In table [5j we give the binding energy, the probability of total nucleon spin being 
equal to two, that for it being equal to one, the charge radius and the root-mean-square 
radius for the ground state and the first excited state. The probability of S = 2 is often 
called .D-state probability, that of S = 1 sometimes P-state probability. Since, however, 
various orbital angular momenta of two can be coupled to a total angular momentum 
of zero, which is non-orthogonal to a pure S"-state, we consider the above definition the 
better one. 

In the first column of table |3j, we give the results of the full calculation. All results for 
the ground state agree reasonably well with those of ref. pi, although the probability of 
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Table 5: Properties of the first two + states in various model spaces 





full space 


coef>0.01 


coef>0.015 


minimal 


ref. i 


^channels 


227 


138 


58 


20 




Etot(MeV) 


-25.910 


-25.276 


-24.557 


-23.203 


-25.86(15) 


P S =2(%) 


10.2 


10.5 


10.1 


9.1 


11.3(1) 


Ps=i(%) 


0.19 


0.18 


0.13 


0.08 


- 


r&{fm) 


1.47 


1.49 


1.51 


1.55 


- 


r r ms(fm) 


1.47 


1.49 


1.51 


1.54 


1.50(1) 


E 0+ {MeV) 


-6.417 


-6.205 


-5.797 


-5.245 




P S =2(%) 


6.86 


6.79 


6.59 


6.34 




Ps=x{%) 


0.06 


0.06 


0.03 


0.02 




r chi fin) 


3.09 


3.10 


3.08 


3.04 






3.02 


3.02 


2.99 


2.97 





finding spin two is a bit smaller, which might be due to slightly different parameters for 
the potential used. The S — 1 configurations contribute about one MeV to the binding 
energy, despite their marginal admixture. Due to the different width parameters used 
for the triton and 3 He fragments, the charge and mass radii of 4 He can in principle be 
different. The differences, however, are too small to show up in table [5[ The energy of 
the first excited state is above the 3 He — n threshold. Allowing for all width parameters 
used for the physical channels, it would occur a few keV above the lowest threshold, the 
3 H — p one. Thus the first excited state would no longer be a reasonable approximation 
to the + resonance, but rather of the lowest threshold. We convinced ourselves that the 
model space chosen yields a ground state energy just 60 keV less bound than in the space 
containing all width parameters. This is about the same amount gained by uncoupling 
the physical channels. Considering, however, how much the model space is increased, the 
gain is just a third of the former. 

If we reduce, however, the number of width parameter to only four, then we lose only 25 
keV for the ground state, but three MeV for the excited state. Therefore we consider the 
model space chosen as optimal. The relatively high energy found for this excited state 
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0.0 2.0 4.0 6.0 

q[fm 1 ] 



Figure 3: Point Coulomb form factor for the first excited state (dashed line) compared 
to the ground state form factor (full line). Also shown is the result for the minimal wave 
function (dotted line). 



might explain the difficulties of shell model calculations |23j to reproduce the position 
given in the compilation [|TJ. 

As one might have expected, the radius of the first excited state is much larger than that 
of the ground state. The probabilities for spin unequal to zero, however, are much lower 



than for the ground state, in contrast to the results found in ref. [23], shading doubts 



on the interpretation of this state as a breathing mode p3fl . Also the form factor for 
the first excited state is quite different from that of the ground state (see fig. |3|). The 
minimum is almost at the same position as for the ground state, but for most momentum 
transfers the form factor of the first excited state is about an order of magnitude below 
the ground-state one. 

14 




Figure 4: Comparison of the point Coulomb form factor of the ground state for the full 
wave function (full line), the 58 channel wave function (dashed line) and the minimal 
wave function (dotted line). 

Since a wave function consisting of more than 200 components, each containing five differ- 
ent width parameters, is not instructive, we tried to find a much simpler one by reducing 
the number of configurations. The most obvious idea of omitting all S = 1 components, 
because of their tiny percentage, was not successful, as we lost well over one MeV in bind- 
ing energy. Therefore we omitted all configurations which had (non-orthogonal) expansion 
coefficients below 0.01. Thus reducing the number of channels by roughly 40 percent, we 
lost only a moderate amount of binding energy, but all the other results remained almost 
unchanged (see second column in table |5|). The change in the form factors is too small 
to be displayed in fig. ^|. Omitting repeatedly all configurations with coefficients not ex- 
ceeding 0.015 we could again cut the number of channels by more than half. The binding 
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energies are now already reduced appreciably, but the other results hardly change (see 
third column in table ||); also the form factors are almost unchanged (see fig. ^). 

This wave function, however, is still too complicated to be used in, say, radiative capture 
calculations. Omitting the last remaining S = 1 channels, which had expansion coeffi- 
cients barely larger than those of the deleted states, we lost again more than an MeV. 
Therefore we changed the strategy of channels to be deleted: We deleted every channel 
if the loss of binding energy was below 70 keV. Employing this procedure we ended up 
with just 20 channels. We call this the minimal wave function, because deleting any ad- 
ditional channel reduces the ground state binding energy by more than one hundred keV 
and, more importantly, the energy of the first excited state by more than 350 keV. The 
results for this wave function are displayed in table [| and figs. [2|, [3] and |j. Here the loss in 
binding energy is large, but all other calculated quantities are still reasonable. Therefore 
we consider this wave function useful as a starting point for more detailed studies, like 



d — d radiative capture [27] . Its structure is given in appendix B. 



5 Partial Wave Analysis 

In the following we will present the results of the resonating group calculation and the 
i?-matrix analysis for all relevant partial waves. Because of the large amount of data we 
restrict the presentation mainly to the elastic scattering phase shifts. Only in cases where 
these seem to disagree, we present some additional information, like Argand plots. The 
R- matrix results are given for an energy up to 7.5 MeV in the center-of-mass of the triton- 
proton channel and correspondingly in the other channels. This energy is well within the 
range of data taken into account, and therefore end-of-data effects should be negligible. 
To give an impression how the calculated phase shifts vary with energy, we display the 
RGM results up to an energy of 12.5 MeV. The scattering calculation, however, is stable 
to the highest energy calculated till now, of 80 MeV. We consider as "physical channels" 
in the RGM calculation the t — p and 3 He — n channels and all [2+2] fragmentations, like 
d — d, d — d and d — d, which can be coupled to the quantum numbers considered. 

Let us start with the partial waves of positive parity. The most elaborate one to calculate 
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Figure 5: Comparison of elastic 1 Sq t — p phase shifts extracted from i?-matrix analysis 
(full dots) and calculated for various model spaces within the RGM (full lines) near the 
3 He — n threshold. The numbers besides the lines denote the dimension of the model 
space. Note that the calculated threshold is a bit lower than the experimental value, (cf. 
table §) 

is for the + -channels, because in this partial wave there is the bound ground state and the 
low-lying first excited state. In fig. |5] we compare the t — p phase shifts extracted from the 
i?-matrix analysis with RGM-calculations for various model spaces in the neighbourhood 
of the 3 He — n threshold. The i?-matrix results show a steep rise till the 3 He — n threshold 
exceeding 100 degrees and a gradual decrease above the threshold. 

For the RGM, we give results for model spaces similar to the ones used for the bound 
state (see table |j). For the simplest model space of 5 physical channels, as described in 
the previous section, we find a repulsive triton-proton phase shift. The 3 He — n threshold 
is hardly detectable. Just around 8 MeV there is a sign of the first excited state (see table 



Adding the 60 distortion channels of [3+l]-fragmentations yields qualitatively the R- 
matrix result: a steep rise up to the 3 He — n threshold, which is a bit lower than in 
reality (compare table ^|), and then a gradual falloff. The possibility of a Wigner cusp 
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Table 6: Resonances extracted from the i?-matrix approach and the RGM calculations. 
In the first column the J n value is given, then the dominant isospin. In the third column 
the dominant structure is indicated, if it is nontrivial, then the energy relative to the t — p 
threshold as given in Q. In the next two columns we give the energies of the latest RGM 
calculation (usually the one with the largest model space) and one, which allowed only 
for the few physical channels as described in the text. All energies are given in MeV. 
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at the threshold cannot be decided from the calculation. From the steep rise one could 
conclude the existence of a nearby + resonance above the threshold. The full calculation, 
which includes 204 distortion channels of various kind (see previous section), reproduces 
the i?-matrix results below the 3 He — n threshold almost perfectly, reaches 100 degrees 
at threshold (which is a bit too low) and then falls off gradually. From this calculation, 
we would position the first excited state in agreement with the i?-matrix just below 
the threshold (see table H). In the previous chapter we showed that the bound-state 
calculation for a model space containing the same structures, but fewer width parameters 
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Figure 6: Comparison of elastic 1 S'o 3 He— n phase shifts extracted from i?-matrix analysis 
(full triangles) and calculated for various model spaces within the RGM (dashed lines) 
above the 3 He — n threshold. The numbers besides the lines denote the dimension of the 
model space used. 

on the relative coordinate of the physical channels, puts the state well above the 3 He — n 
threshold. This means that the state is strongly shifted due to the coupling to the open 
channels, a fact also found in the R- matrix analysis (see ref. [|l|). As was mentioned in 
ref. [||], the 5-matrix pole corresponding to the first excited + state in the i?-matrix 
analysis is quite broad, and appears well above the 3 He — n threshold, at about 3MeV 
relative to the t — p mass. In these respects, it differs considerably from the resonance 
parameters determined by the usual prescription (the multi-level generalisation of which 
we have used), and from the apparent position and width of the resonance in the data. 
However, this is clearly the case of a non-Breit-Wigner resonance, as discussed in section 
3, that appears in the 4 He data primarily as a narrow threshold effect (near the 3 He — n 
threshold), but is in fact a relatively broad structure located higher in excitation energy, 
as it appears in a — a scattering |5j and in most calculations, including the present one. 

This full calculation seems to be the optimal one. Near the deuteron-deuteron threshold, 
however, numerical instabilities develop, because of linear dependencies in the [2+2]- 
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Figure 7: Elastic + phase shifts for all physical two-fragment channels. The data from 
the i?-matrix are denoted by full dots (t—p channels), full triangles ( 3 He — n channels), 
and crosses (d — d channels). The calculated phase shifts are denoted by full lines, dashed 
lines, and dashed-dotted lines, respectively. For fragmentations with more than one chan- 
nel, the quantum numbers are indicated. 

fragmentations. In order to get numerical stability we had to remove most of the [2+2] 
distortion channels. Keeping 133 distortion channels yielded numerically stable results 
for the whole energy range calculated. As can be seen in fig. |^, the deviations from the 
full calculation are only minor ones. 

In fig. ^ we present the elastic 3 He — n scattering phase shifts again from the .R-matrix 
analysis and the various model spaces in the RGM. Starting at the threshold the i?-matrix 
phase shifts are negative and fall to about -80 degrees at 2 MeV. Note that in the following 
all energies are center-of-mass energies above the triton-proton threshold. 
For the physical channels alone, the RGM yields also negative phase shifts, but far from 
the R- matrix ones. Also the calculation including the [3+1] distortion channels is still far 
from the .R-matrix data. The full calculation, however, reproduces the data again almost 
perfectly and the reduction to 138 channels yields only minor modifications at the higher 
energies (see fig. |). 
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Figure 8: Argand plot of the 1 5'o t — p ^-matrix element. The full line represents the S- 
matrix deduced from the .R-matrix in the complex plane, the broken line the one calculated 
within the RGM. The dots on the curves denote steps of one MeV in the center-of-mass 
energy. 

In fig. [f] we display the phase shifts of all two-fragment channels over the whole energy 
range. The beautiful agreement at low energies between the .R-matrix and RGM results 
obviously does not carry on to higher energies. (The RGM calculation is the 138 channel 
one.) 

Let us now consider the d — d phase shifts. The calculated x Sq phases agree reasonably 
well with those from the .R-matrix. The 5 D phases are, however, qualitatively different: 
the .R-matrix phases are negative and reach more than -30 degrees at 7.5 MeV, whereas 
the calculated phase shifts are tiny and positive. We will discuss this point in the final 
section. 

The calculated [3+1] phase shifts start to deviate from the R- matrix results around 2 MeV. 
Whereas the .R-matrix phase shifts for the two channels differ by about 180 degrees, a 
difference which does not modify any experimental datum, the RGM phases for the two 
channels converge to almost the same negative value, however, different from the .R-matrix 
result (see fig. ^). 
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Figure 9: Same as fig. [7|, but for the 1 + phase shifts. 

The Argand plot of the t — p S-matrix elements (see fig. ||) hints at the origin of this 
difference. The diagonal S'-matrix elements become (very) small in both approaches, due 
to the underlying + , T = resonance, discussed above. Since this resonance has almost 
good isospin T = 0, the coupling S'-matrix element between the triton-proton and the 
3 He — n channel reaches almost the unitary limit of unity. Even at 7.5 MeV the coupling 
matrix element exceeds the diagonal ones appreciably in both approaches. Since the 
diagonal S-matrix elements are so small, elastic scattering will not allow to determine 
them more precisely (see discussion below). 

From an analysis of the calculated eigenphase shifts together with the background phases 
for channel radii as given in ref [l], we deduce two additional resonances of almost pure 
d — d structure between 10 and 11 MeV. These resonances do not occur in the recent 
compilation [l[|. The one of 5 D structure appears already in the model space of only 
physical channels (see table |6|). These high lying resonances of pure d — d structure will 
occur in more partial waves and we will discuss them at the end of this section. 

Let us now come to the next partial wave, the 1 + one. In this partial wave the R- 
matrix analysis JIJ finds only one broad resonance about 9 MeV above the t — p threshold. 
Therefore we restricted the RGM calculation to the physical channels only, the 3 Si and 
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Figure 10: Comparison of 2 + t — p and 3 He — n phase shifts. 

3 Di channels for the [3+1] fragmentations and 5 Di for the deuteron-deuteron channel. 
In fig. |9| we compare the elastic phase shifts from the i?-matrix and the RGM for all 
five channels. All .R-matrix phases are negative, with the S-phase shifts reaching -70 
degrees, whereas all D-phases do not even reach -15 degrees. The calculated phase shifts 
are negative too. The S-phases agree nicely with the R- matrix ones and the .D-phases 
are even smaller. From an analysis of the calculated eigenphase shifts together with the 
background phases for radii as given in ref |]J , we again try to deduce resonances. We find 
in the energy range considered solely a single broad T = resonance around 11 MeV (see 
table 0). This is appreciably higher than in the i?-matrix analysis but the resonance 
has the same structure. 

Since we did not allow for distortion channels, which lead to a stronger coupling of channels 
in the interaction region, we anticipate that the energy of this resonance would be reduced 
when distortion channels are included, so that the resonance energies agree more closely. 

The 2 + -partial wave contains the most coupled channels. In the i?-matrix analysis the 1 D 2 
and 3 D 2 [3+1] fragmentation channels are taken into account together with the 5 62, 1 D 2 , 
and 5 D 2 deuteron-deuteron channels. The 5 G 2 d — d channel is neglected. In the RGM 
calculation, we again restrict ourselves to the physical channels and allow for the same 7 
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Figure 11: Comparison of 2 + d — d phase shifts, 
channels. A test calculation [|27| including the 5 G 2 channel, but partially simpler internal 



functions, yields only minor modifications and is therefore not discussed here. In fig. |10 
we compare the R- matrix [3+1] phase shifts with those calculated from the RGM. All 
i?-matrix phases are positive, the 1 D 2 reaching almost 10 degrees, with the t — p phases 
a bit larger than the 3 He — n ones. The 3 D2 phases are quite small, barely reaching 3 
degrees. The l D 2 RGM phase shifts agree nicely with the i?-matrix ones, even so they 
do not reach quite as high. The 3 D 2 phases, however, have essentially the opposite sign. 
It should be noted, however, that a previous i?-matrix analysis using a somewhat smaller 
data basis also yielded negative 3 D 2 phase shifts. The origin of this sign change is not yet 
known. 

In fig. [TTj we present the d — d phase shifts. The 5 S 2 phases fall off strongly and agree 
perfectly for the two approaches. All D-w&ve phase shifts are small. The i?-matrix analysis 
finds positive values at the higher energies, whereas the RGM finds negative values. 

In the recent compilation [0] three 2 + T = levels are given, all essentially of rather pure 
d — d structure, at 7.6, 8.9, and 10.1 MeV above the t — p threshold, having primarily 5 S 2 , 
1 -D 2 , and 5 D 2 components. In the RGM calculation, we cannot easily find the complex 
energy poles of the 5*-matrix, so that the partial widths of the various channels can 
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only be given in case of a Breit-Wigner resonance. We can, however, always diagonalize 
the S'-matrix to get the eigenphases and the eigenvectors at all real energies. From 
these eigenvectors we can determine the isospin value and even its purity of the level 
under consideration and its structure. In the energy range considered, we find also three 
T = resonances, having predominantly d — d structure, at 7.3, 10.7, and 11.0 MeV. The 
structure is the same as found in the R- matrix analysis, but the state at 10.7 MeV has 
also components of X D 2 t — p and 3 He — n of about 30 percent. Since the width of all 
these resonances is some MeV U, the agreement is quite good (see table |]). 

For the 3 + and 4 + partial waves, the .R-matrix analysis finds small negative values for 
all D-wave phase shifts for all fragmentations. Till now there is no full-fledged RGM 
calculation for these partial waves. In order to get an idea what might come out of such 
a calculation, we used all the matrix elements calculated so far that could be coupled to 
3+ or 4+_ p or th e [3+1] fragmentations, the resulting internal wave functions are quite 
good, which can be deduced from the resulting threshold energy, where we loose some 
hundred keV only. The deuteron-deuteron channels, however, are almost unbound relative 
to the four-nucleon threshold in the 3 + case and unbound by more than 3 MeV for the 4 + 
channel. This is also reflected in the calculated phase shifts. The [3+1] phase shifts 
are negative and close to their 3 Di and 3 D 2 counterparts. The 5 D d — d phases are tiny 
and even positive in the 4 + channel. 

Before presenting the negative parity partial waves, let us briefly summarize the results 
for all D-wscves. For the 3 Dj [3+1] phase shifts, the i?-matrix yields negative values for 
J = 1 and J = 3, but positive ones for J = 2. Such a behaviour can be due to a strong 
tensor force. The RGM calculation yields negative values in all cases with barely any 
J-splitting. Also the splitting of the .R-matrix 5 Dj d — d phase shifts follow essentially 
the pattern of a strong tensor force. The corresponding RGM ones are tiny and show 
almost no splitting. 

Let us now come to the negative-parity partial waves. In the compilation || three 0~ 
resonances are given. The .R-matrix analysis takes the 3 Pq for the three physical channels 
into account. In fig. O the .R-matrix phases are compared to the RGM ones. From the 
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Figure 12: Same as fig. [F], but for the phase shifts. 

steep rise of the t — p i?-matrix data around 1.0 and 4.0 MeV, one can easily conclude the 
existence of two resonances. The third one, which is of d — d structure |IJ is not apparent. 
Except where the [3+1] phases vary rapidly with energy, the t — p and 3 He — n phases 
differ essentially by multiples of 180 degrees. Above 4.5 MeV we have added 180 degrees 
to the 3 He — n phase shifts in oder to facilitate the comparison with the RGM results. 

The RGM calculation takes also the physical channels into account and a few distortion 
channels of d — d structure. As for the + channels, the [3+1] results look obviously 
different for the lower energies. At the higher energies they do not quite reach the R- 
matrix phases. The d — d phase shifts agree nicely. Similar to the + partial waves, 
we anticipate that the resonances with their rather pure isospin are the origin of the 
differences. As can be seen from the Argand plots in figs. [TBI and ITJ], this is indeed the 



case. The model space used here corresponds to curve b in figs. [U| and [14]. Even though 
the agreement is not perfect, all the qualitative features agree. Since the i?-matrix and 
RGM curves of the S-matrix pass on different sides of the origin, one of the phases is 
increasing and the other decreasing. From the eigenphase shifts we can deduce also three 
resonances, two of [3+1] structure; a T = one at 1.6 MeV which has a rather large 
T = 1 admixture and a rather pure T = 1 one at 6.6 MeV, and one of d — d structure at 
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Figure 13: Argand plot of the 3 Pq elastic t — p S-matrix elements extracted from the 
i?-matrix (full) and calculated by using the physical channels only (dashed line a), adding 
6 d — d distortion channels (dashed line b) and adding a total of 46 distortion channels 
(dashed line c). The origin is marked as a small circle in the middle of the figure. 

11.1 MeV. These findings agree again nicely with those of ref. |]J (cf. table |6|). 

In order to demonstrate the effect of changing the model space, we display in figs. [13 
and [TJ] the S'-matrix elements for elastic triton-proton and 3 He — n, respectively. For 
the physical channels only, the RGM results looks quite different for the t — p S'-matrix, 
even qualitatively. The standard calculation (curve b) displays already all the qualitative 
features of the i?-matrix results. The model space of dimension 50 (curves c) reproduces 
in both cases the i?-matrix results at low energies perfectly, but crosses the imaginary axis 
below instead of above the real axis (fig. [13]) and vice versa (fig. |j~4]) . From these figures it 
is obvious that small modifications in any of the approaches might change the S-matrix 
a small amount, so that the behaviour near the origin will also qualitatively agree. For 
the large model space, the resonances are all shifted downward to 1.2, 5.8, and 10.9 MeV, 
thus agreeing even better with the compilation [|1|] (see table ||). 
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Figure 14: Argand plot of the 3 Pq elastic 3 He — n S'-matrix elements extracted from the 
R- matrix (full) and calculated by using the physical channels only (dashed line a), adding 
6 d — d distortion channels (dashed line b) and adding a total of 46 distortion channels 
(dashed line c). 



The 1 phase shifts are displayed in fig. [TR The i?-matrix data yield positive 3 Pi and 



negative 1 Pi [3+1] phases and also negative d — d ones. The extracted level structure 
is quite rich ]l|. The RGM calculation uses all physical channels and a few distortion 
channels. The results compare favourably with the P-matrix ones, even though the 3 Pi 
phases are not quite as large. Allowing for a large number of distortion channels, corre- 
sponding to those of the 0~ partial wave, the positive phases grow up to 45 degrees, and 
so come into close agreement. From the eigenphase shifts we deduce also four resonances 
(see table |B|), which agree nicely in structure and position with the P-matrix ones. It 
should be noted that the 3 Pi [3+1] resonances are strongly isospin mixed. 

In the 2~ partial wave the P-matrix allows for P- and F-waves in all configurations. In 



fig. [Tg the P-wave phase shifts are displayed. The [3+1] phases reach about 90 degrees at 
rather low energies, which indicate the existence of two lying resonances (see ref. [l| and 
table |B|). The d — d phases are negative as in the other two negative-parity partial waves. 
The RGM calculation yields somewhat different results: the positive phase shifts reach 



28 



-50.0 



-100.0 1 1 1 1 ' — 

0.0 5.0 10.0 

E c JMeV] 



Figure 15: Same as fig. [7], but for the 1 phase shifts. 



Hililliliilil 

i' , 
















\ 






* — «. — 



0.0 5.0 10.0 

E c JMeV] 



Figure 16: Same as fig. 0, but for the 2 phase shifts. 

only half the i?-matrix value, whereas the negative d — d phases are about a factor two 
more negative. The extracted resonance positions are still quite reasonable (see table 

Since the i?-matrix analysis allows also for F-waves, there are also results for the 3~ and 
4r partial waves. The l F% phase shifts are negative and reach -4 degrees at the highest 
energy. The 3 Fj phase shifts are all positive or tiny, i.e. in magnitude below 0.2 degrees. 
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The [3+1] 3 F 3 phase shifts reach up to 7.0 degrees and the d — d ^F^ one up to 4.8 degrees, 
all others are tiny. Resonant structures do not appear. No RGM results exist till now for 
these partial waves. 

For the negative parities the i?-matrix and RGM phase shifts could be interpreted as due 
to an effective strong tensor force for the [3+1] fragmentations. Whereas the .R-matrix 
results for the d — d F-waves point more into the direction of a stronger spin-orbit force. 

The high lying resonances of d — d structure, which occur in almost all partial waves 
deserve some further discussion. They are not related to any circling of the corresponding 
S'-matrix element as in the + and 0~ partial waves (see figs. |^, [13], and [13). They are due 
to the prescription given in ref. [l], where the d — d channel has a channel radius of 7.0 fm. 
This large channel radius leads to rather large background phase shifts, so that around 11 
MeV the differences between the eigenphases of the S"-matrix and the background phase 
pass through ninety degrees, the criterion we chose for defining the energetic position of 
a resonance. So it might be that the position of all of these high energy d — d resonances 
is strongly influenced by the large value of the channel radius. Unfortunately within the 
RGM we are not aware of any better criterion. This might be the reason, why we find in 
the RGM calculation also in the + and 1 + partial waves resonances around that energy, 
which have no counterpart in the i?-matrix analysis (cf. table |6]). 

6 Comparison with data 

In this chapter we compare the various calculations with a wide variety of experimental 
data. We present figures for all possible elastic scatterings and reactions. Out of the 
large amount of data, we have selected those energies for which several observables are 
measured, if possible by different groups. Most of the data presented here are in the 
data set used to determine the i?-matrix parameters. We will point out those that are 
not included in the set. We present results for the i?-matrix analysis together with 
RGM calculations. In order to demonstrate the effects of enlarging the model spaces in 
the calculation we display curves for the largest model space used, denoted by "full", 
and for the negative parities for the physical model space only, denoted by "small". 
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To demonstrate how the calculation using a realistic nucleon-nucleon interaction differs 
from previous calculations, we display also the results from a previous calculation ||], 
denoted by "semi- realistic" . Since the i?-matrix analysis uses most of the data as input 
values, we consider the S'-matrix elements determined in this analysis as experimentally 
determined, despite the fact that the analyis has not yet fully converged. We emphasize 
the differences of the S'-matrix elements originating from the .R-matrix analysis and the 
full RGM calculation and present the results when we change the value of a specific matrix 
element from its RGM value to its i?-matrix value in order to demonstrate that it is just 
this single matrix element that does not allow to reproduce the data. At the end of this 
chapter we discuss the various differences and their possible origin. It should be noted 
that the partial wave analysis of the RGM calculation and the .R-matrix differ. The R- 
matrix analysis contains additionally the positive-parity contributions for J 71 = 3 + and 
4 + and all the negative-parity F-wave contributions. We will point out in the following 
in which reactions these differences play a substantial role. The meaning of the lines in 
all figures is the same as given in fig. |18| . 

The various reactions are presented in the order of their thresholds, starting with triton- 
proton elastic scattering. Here differential cross section and analysing power measure- 



ments for the proton and also for the triton exist around 4 MeV proton energy. In fig. [L7 
we compare the differential cross section data with the various calculations. 
The .R-matrix analysis is performed in 250-keV steps of the center-of-mass energy starting 
from the triton-proton threshold. Since there are no rapid changes with energy of the 
observables presented, we use the results of the calculations for the nearest energy to that 
of the data. The error introduced by this procedure is usually well within the experimental 
error bars. The RGM calculations are also in 250 keV steps of center-of-mass energy, but 



the calculated thresholds are slightly different (see table |2|). From fig. [17] we see that 
the .R-matrix analysis reproduces the data quite well. This demonstrates how good the 
fit to the data is in general. The RGM calculation, on the other hand, reproduces the 
data reasonably well, taking into account that the calculations contain no adjustable 
parameters. There are, however, marked deviations from the data. The full calculation 
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Figure 17: Differential elastic proton-triton cross section calculated at 3.0 MeV E cm . 
The data are for 4.15 MeV protons from ref. [^8| (dots) and from ref. |29| (triangles). 
The full line represents the .R-matrix analysis, the dashed one the full calculation, the 
dot-dashed one the small calculation, and the dotted one the semi-realistic calculation. 
The errors in the data are of the size of the symbols. 

using the Bonn potential describes the data much better than the small calculation, but 
the deviations for forward angles are appreciable. 

The calculation using the semi-realistic potential describes the forward data very well, 
but deviates strongly for larger angles. Since we know the S'-matrix elements from the 
calculations and the .R-matrix analysis themselves, we can try to trace the differences 
to specific ^-matrix elements. The full calculation differs from the .R-matrix analysis 
in the 3 P 2 matrix element. Note that we represent all S'-matrix elements in the form 
Ski = Vki£ 2l5kl with phase shifts Ski and channel couplings 77^. In the full calculation, 
77 = 0.86 instead of 0.68 for the .R-matrix analysis, and the phase shift of 49 degrees 
misses by 20 degrees the R- matrix result; thus the coupling to the 3 He — n channels is 
much weaker. When we use the .R-matrix 3 P 2 matrix element instead of the calculated 
one, the data are very well reproduced. The results lie on top of the .R-matrix results, 
therefore we do not display them in fig. [17|. The semi-realistic calculation yields a much 
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Figure 18: Proton analysing power of the reaction 3 H(p, p) 3 H calculated for 3.0 MeV E cm . 
The data are for 4.15 MeV protons from ref. The full line represents the i?-matrix 

analysis, the dashed one the full calculation, the dot-dashed one the small calculation, and 
the dotted one the semi-realistic calculation. The open circles denote the full calculation 
with the 3 P 2 matrix element replaced by the corresponding i?-matrix one. 

better 3 P 2 matrix element. Here the deviations are mainly due to too strong a coupling 
for the 3 Pi matrix element. 

Proton analysing power measurements at the same energy as the cross section and nearby 



triton analysing power data are displayed in figs. [18| and [19], respectively. The i?-matrix 
analysis again does an excellent job, whereas the RGM calculation reproduces the data 
only qualitatively (see figs. |l^ and [19]). This is again due to the differences in the 3 P 2 
matrix elements (see open circles in figs. [18] and [19]) , which reproduce the data much better. 
Note, however, that all calculations yield an analysing power for the proton and the triton 
that is markedly different in the maximum. This difference is due to the coupling matrix 
element 1 P\ — ► 3 Pi, which is different from zero due to the spin-orbit component of the 
potentials used. Contrary to the nucleon-nucleon system, such singlet-triplet transitions 
are not forbidden by the Pauli principle. 

The differential cross section for the 3 H(p, n) 3 He reaction at almost the same energy is 
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Figure 19: Triton analysing power of the reaction 3 H(p,p) 3 H calculated for 3.0 MeV E, 
The data are for 4.28 MeV protons from ref. pOf. The labeling is the same as in fig. [18 
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Figure 20: Differential cross section of the reaction 3 H(p, n) 3 He calculated for E cm = 3.0 
MeV. The data are for 4.101 MeV protons from ref. |3T| . The labeling is as in fig. [IB]. 

displayed in fig. |2(]. Whereas the i?-matrix reproduces the data quite well, the RGM 
calculation yields only qualitative agreement with the data, with the increased model 
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Figure 21: Differential cross section calculated for 8MeV neutrons off 3 He. The data are 
for 7.9 MeV neutrons from ref. P31 (dots) and from ref. [B2J (triangles) for 8 MeV. The 



labeling is as in fig. [TJ. 



space yielding again better agreement. The differences can again be traced to essentially 
a single matrix element, the already known 3 P 2 one. For the Bonn potential it reaches 
only two thirds of the magnitude of the P-matrix one, and its relative phase to the 1 Sq 
matrix element is only 60 degrees instead of 80. This just demonstrates the unitarity of 
the S'-matrix, i.e. the missing coupling in the elastic triton-proton channel now shows up 
in the missing strength going from the proton to the neutron channel. All other large 
matrix elements agree within a few percent in modulus and phase. Changing again in 
the full calculation the 3 P 2 matrix element to its P-matrix value reproduces the data 
almost perfectly (cf. fig. |2(|). Except in the minimum, the modified RGM agrees with the 
P-matrix result. For the semi-realistic potential, the 3 P 2 matrix element overshoots the 
P- matrix one by 25 percent and the 3 Pi one by even 50 percent. 

Elastic neutron scattering cross section and neutron analysing powers around 8 MeV 
neutron energy are displayed in fig. [ZT] and fig. respectively. The cross section data 
are well reproduced not only by the P-matrix analysis but also by the full and semi- 
realistic calculations. Changing in the full calculation again the 3 P 2 matrix element to its 
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Figure 22: Neutron analysing power for the scattering of 8 MeV neutrons off 3 He. Data 



are from ref. (triangles) and from ref. |35|] (dots). The labeling is as in fig. IS. 
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Figure 23: Differential cross section for the reaction 2 H(d, p) 3 H calculated for E cm = 2.11 
MeV. Data are for 4.0 MeV deuterons from ref. [BBf (triangles) and from ref. [37| (dots). 



The labeling is as in fig. IS 



i?-matrix value brings the results on top of the semi-realistic cross section results. 
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Figure 24: The analysing power iT u for the reaction 2 H(<i,p) 3 H calculated for E cm = 2.11 
MeV. Data are for 4.0 MeV deuterons from ref. [BBf.The labeling is the same as in fig. ITB 



The neutron analysing power, however, is well reproduced by the i?-matrix, but only 
qualitatively by the RGM-calculations. The height at the maximum analysing power 
around 115 degrees is controlled by the phase of the 3 P2 matrix element, since the modulus 
of the matrix elements turns out to be quite similar for all calculations. In the full 



calculation the phase shift is about 25 degrees below the .R-matrix analysis (see fig. [T6p . 
Increasing the phase of the RGM matrix element by this amount, as was done above for 
the cross section, reproduces also the analysing power data perfectly. Also the minimum 
around 90 degrees is well reproduced for this choice. It should be mentioned, however, that 
the coupling matrix element x Pi — > 3 Pi is only about half the P-matrix one, thus leading 
to a much smaller difference between the neutron and 3 He polarisations (compare with 
the discussion of the charge conjugate scattering 3 H(p, p) 3 H above). There are, however, 
no 3 He polarisation data available, hence one cannot decide about which calculation is 
correct. All the above elastic neutron scattering data are not included in the i?-matrix 
fit, hence, we consider the calculations to be predictions. 

Let us now come to the deuteron induced reactions. In fig. |23] the 2 H(d, p) 3 H differential 
cross section is displayed. The data are well reproduced by all calculations. Here for the 
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first time the small calculation reproduces the data even better than the full calculation. 
The semi- realistic calculation overshoots the data a bit. At the forward (and backward) 
angles, the F-waves, which are missing in the RGM calculations, contribute about 50 
percent to the cross section in the i?-matrix analysis and indicate that even G-waves may 
be necessary. 



In figs. 0, |25], |26|, and \27\ the deuteron analysing powers are displayed. Whereas the 
i?-matrix analysis describes the data nicely, the RGM calculations fail completely. In the 
i?-matrix analysis the dominant matrix elements are the 3 Pi, 1 D 2 , and the 5 S 2 — > 3 D 2 
matrix elements, but also quite a number of others must not be neglected, like the 1 Sq 
and the 5 S 2 — > 1 D 2 . The coupling matrix element 5 S 2 — > 3 D 2 determines essentially the 
form of the angular distribution of all polarisation observables. The F-matrix analysis 
finds this matrix element to half of the leading 3 Pi one and the relative phase shift to 
be +60 degrees. The full RGM calculation yields, however, only a quarter of the leading 
matrix element and the relative phase shift of -33 degrees, i.e. the opposite sign. The 1 So 
and 3 P 2 overshoot the i?-matrix result by a factor 2 to 3. The other RGM calculations 
give similar results. Changing the RGM coupling matrix element 5 S 2 — > 3 D 2 to the R- 
matrix value yields reasonable agreement for the vector polarisation data (cf. fig. |24|) and 
improves the angular form of the tensor analysing powers (cf. figs. |26], and ^). For 
the forward angles the calculation comes now much closer to the data, but the shape of 
the angular distribution does not match. The reason for this discrepancy can be found 
in the interference terms of the large and the smaller matrix elements. The P-matrix 
analysis indicates that for a detailed description of these data, the F-wave contributions 
are necessary. The modified cross section turns out to be slightly too large. 

For the 2 H(<i, n) 3 He reaction, data exist for the same energy as for the charge conjugate 



2 H(<i,p) 3 H reaction. In fig. ^8| we compare the differential cross section data with the 
various calculations. The i?-matrix analysis again does a perfect job. Also the full RGM 
calculation reproduces the data nicely, whereas the semi-realistic one overshoots the data 
appreciably. The dominant matrix elements are again 3 P 1; 1 D 2 , and 5 S 2 — > 3 D 2 in the 
i?-matrix analysis, but also the 1 So, 3 P 2 , and 5 S 2 — > 1 D 2 being of equal magnitude, must 
not be neglected. The full RGM calculation yields the two dominant matrix elements 
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Figure 25: Same as fig. E3, but for the analysing power 
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Figure 26: Same as fig. |24|, but for the analysing power T 2 \. 



similarly, but the coupling matrix elements 5 S 2 — > 3 D 2 and also 5 S < 2 — > 1 D 2 are only half 
of the .R-matrix value. On the other hand, the 1 Sq and 3 P 2 ones are doubled, as is the 
case of the charge conjugate reaction. 

The Ohio-State group [R3] has measured all deuteron polarisation observables also for 4.0 



MeV deuterons. In figs. 2S and 3C, we present only the vector analysing power and A z 
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Figure 27: Same as fig. 24, but for the analysing power T 22 - 
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Figure 28: Differential cross section for the reaction 2 H(d, n) 3 He calculated for E cm = 2.11 



MeV. Data are for 4.0 MeV deuterons from ref. J37J. The labeling is as in fig. [18 



data, because the overall behaviour is quite similar to the charge conjugate reaction. 
As in the charge conjugate reaction 2 H(<i,p) 3 H, only the i?-matrix analysis reproduces 
the data, whereas all RGM calculations fail completely. To find such a good agreement 
with data, the 3 F 3 and 3 i<4 matrix elements in the .R-matrix are vital. Again the 5 S 2 — > 



40 



0.4 
0.3 
0.2 
0.1 



2 H(d,n) 3 He 




180 



cm. 



Figure 29: The deuteron vector analysing power for the reaction 2 H(d, n) 3 He calculated 
for E cm = 2.11 MeV. Data are for 4.0 MeV deuterons from ref. f38|.The labeling is the 



same as in fig. 17 
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Figure 30: Same as fig. |2~9], but for the analysing power A zz . 



3 D2 matrix element is dominating the angular form of all analysing powers. Using in 
the full RGM calculation the corresponding 5 S , 2 — > 3 D 2 matrix element from the R- 
matrix does not reproduce the data much better. Here it is also neccessary to modify the 
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Figure 31: Differential cross section for elastic deuteron-deuteron scattering calculated 
for E cm = 3.11 MeV. Data are for 6.0 MeV deuterons from ref. [[39]. The labeling is as in 



fiff. 17 



corresponding 5 S*2 — > 1 D2 matrix element. Then the cross section is not changed at all, 
the vector analysing power is reproduced reasonably well, but for the description of the 



tensor polarisations, not too much is gained (cf. figs. |29| and |30[). Note that in the RGM 
no F- waves are included. 

Let us now discuss the last two-fragment process, elastic deuteron-deuteron scattering. 
Because of the identical bosons in the entrance and exit channel, the vector polarisation 
and T21 are antisymmetric and all other observables are symmetric about 90 degrees. Since 
the data are usually converted into the forward hemisphere we display this part only in the 
following figures. In fig. |3I| we present rather old data [^] for the differential cross section 
together with the calculations. The i?-matrix analysis reproduces the data very well. 
The two realistic RGM calculations lying on top of each other also reproduce the data 
reasonably well, the semi-realistic calculation, however, overshoots the data appreciably. 



The analysing powers have been measured by the Zurich group 



In figs. 32, B3L 34 



and |3S| we compare the data with various calculations. All polarisations are quite small, 
of the order of one or two percent. Note the large scale on the figures. Except for the 
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Figure 33: Same as fig. |32] but for the tensor analysing power T 20 . 



tensor polarisation T22, where all calculations fail, the i?-matrix analysis reproduces the 
data well. Considering the small values, also the full calculation reproduces the data 



43 




-0.05 



30 



60 



90 







cm. 



Figure 34: Same as fig. B3, but for the tensor analysing power T 2 \. 



reasonably well. Since the polarisations are so small, it is useless to search for the matrix 
element (s), which dominate the structure of the polarisations. It has to be mentioned, 
however, that the inclusion of these tiny polarisations into the i?-matrix analysis changed 
the sign in the 3 D 2 triton-proton and neutron- 3 He elastic matrix elements (see the above 



discussion following fig. |10|) . 

Considering all the figures displaying a comparison of data and calculations together, it 
is obvious that the parameter-free calculations reproduce the differential cross sections 
rather well, whereas the polarisations are sometimes grossly missed. The differences in 
the description of the data between the i?-matrix analysis and the RGM-calculations are 
concentrated in a few S'-matrix elements. In order to reproduce the data the calculated 
3 P 2 matrix elements for all channels and the 5 S 2 — > 3 D 2 and the 5 S 2 — > 1 D 2 matrix 
element in the deuteron-deuteron fusion reaction have to be changed. In the 3 P 2 one, 
the coupling to the F-waves cannot be totally neglected. So in both cases the tensor 
force is crucial for the determination of these matrix elements, because the coupling is 
to channels with a different orbital angular momentum. As discussed in the previous 
chapter, the i?-matrix analysis indicates a rather strong tensor force. It might therefore 
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Figure 35: Same as fig. 32, but for the tensor analysing power T22. 



be that the Bonn potential is the origin for the discrepancies of the RGM calculation with 
data. That the tensor component in the Bonn potential is rather weak has been pointed 



out many times, e.g. [41 . 



7 Conclusion and Outlook 

We performed microscopic multi-channel calculations employing an r-space version of 
the realistic Bonn potential in the framework of the resonating group model. This is a 
parameter-free calculation in a rather large model space, chosen by physical arguments 
and restricted by computational limitations. For the ground state it is essentially a 
variational calculation, with the functional form of the wave function dictated by the 
following scattering calculation. 

For the ground state energy we got almost the same results as were obtained using varia- 
tional wave functions || for the same potential. The extracted point Coulomb form factor 
agrees nicely with that of a variational calculation p3] using another realistic nucleon- 



nucleon potential. Restricting the model space to those parts appropriate for bound state 
wave functions, we could determine properties of the first excited state and also the point 
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Coulomb form factor. From the properties we concluded that this resonance is not a 
breathing mode and gave arguments why the energy of this resonance in shell model cal- 
culations differs from that determined from scattering experiments. Also the form factor 
deviates strongly from that of the ground state. 

From the scattering calculations we extracted phase shifts and demonstrated that the 
parameter-free calculation yields not only qualitatively, but in most cases even quanti- 
tatively, the same answers as the .R-matrix analysis using only data. In some cases we 
showed that the apparent differences of the two methods are caused by only minor vari- 
ations of the 5*-matrix elements near the origin. The other deviations appear in channels 
that can be strongly influenced by the tensor-force terms in the nucleon-nucleon interac- 
tion. Since it is well known that the potential we used has a rather weak tensor component 



Hfl , we blame the potential for those differences. Extracting resonance positions in an 
analogous way to the recent compilation [[[]] yielded close agreement in the sequence, quan- 
tum numbers, and structure of the resonances. Near the upper limit of the energy interval 
considered, the difference between eigenphases of mainly deuteron-deuteron structure and 
background phase shifts passed through ninety degrees, the accepted criterion for the po- 
sition of a resonance. In these cases, however, the interpretation as a resonance seems to 
be doubtful. 

Finally we compared the various calculations with data. All differential cross sections 
were well reproduced, especially if we used the largest model space. The polarisation 
observables, however, revealed the shortcomings mentioned above. Only after modifying 
in most cases a single S'-matrix element, the agreement between calculation and data 
became satisfactory. For the d — d fusion reaction, however, F-wave contributions that 
could not be taken into account and the interference effects with smaller matrix elements 
prevented a reasonable agreement. 

Taking all this information together, the parameter-free microscopic calculations using a 
realistic nucleon-nucleon force are very promising in explaining and even predicting data. 
The shortcomings discussed above indicate the direction of further studies. Increasing 
the model space brings calculations and data in closer agreement (cf. figs. |T7| to BH), but 
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the essential deviations between i?-matrix analysis and RGM calculations are untouched. 
Since we suspect the nucleon-nucleon potential is responsible for these deviations, it is 
necessary to use another potential with a stronger tensor component. Since the RGM 



requires the potential to be given in r-space, we consider the Argonne AV14 [42] the 
proper choice, especially because there exist many calculations of ground state properties 
for light nuclei with this potential. Furthermore, results for this two-body potential 
supplemented by various three-body potentials are also available. Work in this direction 
is under way. 
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A Appendix 1: Wave functions of the fragments 

The general structure of the wave functions is given in eqs. (1) and (2) (for further details 
see j8], 0). For the deuteron we have used the most simple wave function. Since the spin 
part is trivial, we give just the orbital part. 

= E c n e~^ 2 ^~^ 2 y (r x - r 2 ) + E c n e~^ 2 ^~ r ^ y 2 (r 1 - r 2 ) (14) 

n=l n=4 

Here y denote solid spherical harmonics []57| . The width parameters (3 n and the coeffi- 



cients c n for the deuteron and the singlet deuteron d are given in table [7|. 
The spatial part of the three-nucleon wave functions is of the form 

^n= E c^ L y h ( Sl ) y h (s 2 ) <P S (15) 

h,l2,L,n,m 



47 



where cf) S denotes the appropriate spin functions of each term in the sum. The Jacobian 
coordinates are given by 

1 . 1 . . [2 

Si = -^=(n - r 2 ) and s 2 = -^(ri + r 2 ) - y - r 3 (16) 

The width parameters (3 n and j3 m for triton and 3 He are given in table The first 
three parameters belong to the structures containing no D-waves and the last three width 
parameters belong to those containing at least one D-w&ve. 



B Appendix 2: Minimal 4 He Wave Function 

Since the full bound state wave function of the 4 He ground state contains 227 different 
configurations, it is too complicated to be given in detail. To give an impression about the 
dominant structures of the ground state wave function, we will describe in some detail the 
minimal wave function of table [|. It consists of 20 configurations, where 11 originate from 
the physical channels, which have been resolved into its constituents, and 9 are former 
distortion channels. All of them have orbital angular momentum zero on the relative 
coordinate between the fragments, but as can be seen in table ||, all spin values are present. 
For the relative motion we choose 5 width parameters, used also for all distortion channels 
irrespective of the fragmentation, with the values: 2.433, 0.8214, 0.2987, 0.1185, 0.05001 
in units of fm -2 . For the internal width parameters of the fragments the parameters given 
in tables [7| and || are used. 

The dominant S = contribution contains 11 configurations: In the triton-proton frag- 
mentation, the three containing no internal angular momentum from the physical channel 
and also the combination with the proton-neutron pair inside the triton coupled to spin=l, 



Table 7: Deuteron and singlet deuteron parameters 



n 


1 


2 


3 


4 


5 


Pn 


5.8359 


0.5201 


0.0684 


1.625 


0.4339 


Cd 


6.239420 


-6.659597 


-2.370518 


-6.685710 


-0.7940629 


Cd 


-4.491821 


3.059319 


3.018020 
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Table 8: Triton and 3 He width parameters 



k 


1 


2 


3 


4 


5 


6 




4.972 


0.8659 


0.1531 


8.665 


1.8056 


0.5414 


ym 


3.293 


0.3972 


0.07513 


1.626 


0.4745 


0.09417 


8 3Re 


4.8887 


0.8286 


0.1493 


8.354 


1.939 


0.61723 




3.616 


0.3822 


0.06763 


1.546 


0.4668 


0.09258 



but with the coefficients for the various internal width parameters taken for the third ex- 
ited state. The 3 He — n fragmentation contains the analogous three contributions from 
the physical channel. From the distortion channels we also find the vector with the in- 
ternal proton-neutron pair of S = 1, but here the second and third excited states appear 
with large expansion coefficients, indicating that together with some other state they are 
numerically almost linearly dependent. For the deuteron-deuteron fragmentation none of 
the physical components contribute. The S'-wave part of the deuteron, together with the 
S-wave part of the first and second excited states of the deuteron, form the remaining two 
channels. The overlap of the components originating from the physical channels with the 
total wave function is in excess of 90 percent. The overlap of the configurations containing 
excited fragment wave functions is only about 84-85 percent, whereas the one containing 
the first excited state of the deuteron has less than 10 percent. 

The next important S — 2 contribution contains 8 configurations: For the triton-proton 
fragmentation only the two configurations with one internal angular momentum being 
two, the other zero, originating from the physical channel, contribute. For the 3 He — n 
fragmentation, the configuration with P-waves on both internal coordinates and that 
with a D-wave between the internal neutron-proton pair contribute from the physical 
channel. From the distortion channels the two configurations with one internal D-wave, 
with the coefficients from the first excited state contribute. For the deuteron-deuteron 
fragmentation, the S'-wave part of the deuteron and the D-wave of the deuteron and the 
first excited state of the deuteron form the two contributing channels. The overlap of the 
various components with the total wave function varies between 3 percent for the d — d 
distortion channel and 27 percent for the 3 He D-wave physical channel. 
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The 5 = 1 component is formed from the D-wave part of the deuteron and the D-wave 
of the first excited state, coupled to spin equal one. 
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